Connecting the discrete and continuous-time quantum walks 
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Recently, quantized versions of random walks have been explored as effective elements for quantum 
algorithms. In the simplest case of one dimension, the theory has remained divided into the discrete- 
time quantum walk and the continuous-time quantum walk. Though the properties of these two 
walks have shown similarities, it has remained an open problem to find the exact relation between 
the two. The precise connection of these two processes, both quantally and classically, is presented. 
Extension to higher dimensions is also discussed. 
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Continuous-time quantum walks (CTQW) were intro- 
duced by Farhi and Guttman [3| as generalizations of 
diffusion-type differential equations, in which probabil- 
ity is replaced by a complex amplitude and Markovian 
dynamics is replaced by unitary dynamics. Their mo- 
tivation was to explore whether, in a given framework, 
coherent quantum processes could show dramatic differ- 
ences from classical random walks. This turned out to 
be the case, for Childs et al. demonstrated a graph prob- 
lem that could be solved using the CTQW exponentially 
faster than not just classical random walks but all classi- 
cal methods In one dimension, the CTQW is simply 
the finite-difference Schrodinger equation [2j 

idtip{n, t) = -7[-0(n -I- 1, i) - 2V'(n, t) + - 1, (1) 

where '0('^7 1) is a complex amplitude at the (continuous) 
time t and (discrete) lattice position n. 

The discrete-time quantum walk (DTQW), introduced 
by Aharonov et al. 4] and independently by Meyer [^, 
is a discrete unitary mapping such as 



iiRln, T 4- 1) = cos O^iRin - 1, r) 
V'L(n, r + 1) = cosOiJjLin + 1,t) ^ 



' i sinOipL{n — 1, t) 
i sin 9ipji{n -\- 1, r), 

(2) 

where ipR^n^r) and ^pL{n,T) are complex amplitudes at 
the (discrete) time r and (discrete) lattice position n, 
and the labels R and L indicate an additional degree of 
freedom, often taken as the state of a coin which tells the 
walker (located at position n) which way to step. This 
discrete dynamics has a rich mathematical structure that 
is quite foreign to the CTQW, and has been the subject 
of extensive theoretical investigation. In particular, there 
has been significant extension of the DTQW to include 
decoherence Q, quantum chaotic and quasiperiodic 
effects 1^. In addition, certain continuum limits have 
been used to connect the DTQW to more familiar wave- 
like propagation (f^ [ii| . Despite this large body of 
work, the relation of these two quantum walks remains 
an open problem. 

This problem is truly fundamental for quantum com- 
putation, for at least two reasons. First, it is quite unnat- 



ural to have two distinct ways to quantize classical dif- 
fusion. Determining whether quantum mechanics speeds 
up a classical process is difficult enough, but even moreso 
if there is no unique quantization. For the processes con- 
sidered here, the coin degree of freedom appears unneces- 
sary, and indeed there is a perfectly reasonable discrete- 
time quantum process that can be implemented without 
a coin 12] — this will be discussed below. Second, the 
spreading properties of the two quantum walks are quite 
similar From an initially localized state, both evo- 
lutions generate a probability distribution that is nearly 
constant save for two peaks at ztct (here c = 27 for 
and c — cos 9 for (O), decaying to zero thereafter (see 
Figure 1). Both have standard deviations of position that 
grow linearly in time, quadratically faster than classical 
diffusion. These similarities suggest that, besides the fact 
that both are unitary quantum processes, there should be 
some underlying connection between the two walks. Nev- 
ertheless, the precise relationship has remained elusive. 
In particular, no one has demonstrated how to get ^ 
from (PI by some limiting process. 

Aside from certain approximations verified numerically 
[nl Hsf . the closest previous connection of these two 
walks is the weak limit theorems for the probability den- 
sity due to Konno and others 0|. Specifically, letting 
n ^ a;, where x is considered a continuous variable, the 
long-time limit of the CTQW probability density is 



PcTQw{x,t) 



I 



(3) 



where — 27< < a; < 27<, while the long-time limit of the 
DTQW probabihty density is 



sin 9 



7r(l — x'^T^'^)^/(coa'9Ty^'^^ll? 



(4) 



with — cos 9t < X < cos 6t. Comparing these two ex- 
pressions, one might be led to consider the limit r — > 00, 
6 7r/2, such that cos 9t 2^t. This certainly maps Q 
to 0, but what about and (O? Here I will show that 
this hmit does indeed map the DTQW to the CTQW, has 
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(a) DTQW (b) Continuous-time limit of (c) CTQW 
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FIG. 1: Evolution of (a) the DTQW with cos 61 = 1/4, (b) the 6 tt/2 continuous-time limit of the DTQW, with 7 = 1/8, 
and (c) the CTQW with 7 = 1/8. The probability density p{n,t) — '^{n,t)^'^{n,t) is shown, where the initial conditions for 
(a) and (b) are given by l|16|l. while for (c) ip{n, 0) = 5n,o- 



a direct parallel with the relevant classical random walks, 
and can be extended to higher-dimensional walks. 

The most well known limit of the DTQW with 
6* ^ 0, was introduced by Feynman in 1946 to construct 
apath integral for the propagator of the Dirac equation 
|l5|. In Feynman's picture, a particle zig-zags at the 
speed-of-light across a space-time lattice, flipping its chi- 
rality from left to right with an infinitesimal probability 
each time-step. The resulting dynamics, in the contin- 
uum limit [231, is the Dirac equation, with the flipping 
rate determined by the mass of the particle. 

The limit considered here is, at first sight, quite puz- 
zhng. With 9 7r/2, the DTQW describes a particle 
flipping its chirality with nearly unit probability each 
time-step. Such a particle should not move at all, and 
indeed the maximum group velocity is cos 9 [l]| . which 
goes to zero in the same limit. However, one must also go 
from discrete to continuous time — only by taking these 
two limits together does the DTQW become the CTQW. 

First, it is convenient to work in momentum space by 
introducing the Fourier transform 

(5) 

Then, the DTQW after r steps is the unitary mapping 

where U (t) — If^ with the single-step propagation ma- 
trix 

rr -tka, -tea^ _ { e-^^cos9 -16-''= sin 61 \ , , 
\-ie^''sm.9 e^^cos9 J' ^' 



and I have introduced the Pauli matrices <Jx,<^y,<^z- 

The key step is to set 9 = tt /2 — 5, where (5^1. Then, 
using e"**'^^ = —iaxe^^"^^, I compute U'^: 

— j"j2g— ifeo-j gi(5(T^ gifco-^ giicr^ 

= (— j)^ exp[i(5(cr^ COS 2fc + ay sin 2fc)]e"''^=" 
= (-i)2exp[i(5(cr:r(l-fcos2fc) + <jy sin2fc)] -f 0((5^) 
= {—if' exp[i(52 cos k{ax cos k + ay sin k)] + 0{6^), 

(8) 

where I have used properties of the Pauli matrices in the 
second and third lines, the Bakcr-Campbell-Hausdorff 
theorem in the fourth, and trigonometric identities in 
the last. Now, applying U'^ for t/2 times, and taking the 
limit (5 ^ 0, T(5 -> 27t, I find 

U{t) = (— i)'^ exp[i27t cos /c((J2: cos fc -|- (Tj, sin fc)] , , 
= exp(-i$) exp(-iiJi), ^ ' 

where I have defined $ — tit/ 2 and 

if = —27 cos fc(crx cos fc -I- cr sin fc). (10) 

Using the result lO in © shows that aside from the 
unimportant overall phase <&, this limit of the DTQW 
describes continuous time evolution with a Hamiltonian 
given by (|10|l . The corresponding Schrodinger equation 
idt'i' = H'^, found by the Fourier transform (O, is 

idttpR{n,t) = -7[V'L(n,t) + V^L(n - 2,i)] 
idti^L{n,t) ^ --f[ijR{n,t)+i:R{n + 2,t)], ^ ' 

where here and in the following the amplitudes have been 
trivially redefined, e.g. ipii(n,T) 'iJjR{n,t). 

The reduction to the CTQW is nearly complete. To 
make the connection fully transparent, first observe that 



3 



a general solution of can be split into two terms 

= e*^''**+(n,t) + e-''2'^**_(n,t) (12) 

where 



ipL{n,t) + V'fl(«+ 1,*) 



and 



By direct substitution one finds 
idt'i/±{n,t) 

= T7[*±('^ + 1, ^) - 2*±(n, t) + ^±{n- 1, t)]. 

(15) 

Thus, the continuous-time limit with — > 7r/2 of the 
DTQW is really equivalent to two copies of the CTQW. 
As a consequence, there still remains interference as- 
sociated with the combination of the two components 
'^±{n,t) in H12|l . which oscillates in time. This can be 
eliminated by requiring ipn^rijt) = ipLin — 1,^), so that 
^'_(n,i) = 24]; an equivalent simplification was pre- 
viously found by projecting the initial state onto the 
"positive-energy" subspace of U 

The simplest initial condition with 5'_(n,t) = and 
localized symmetrically about n = takes the form 



V'i?(",o) 



1 / Snfl + Sn,l 

2 1 Sn-1 + Snfl 



(16) 



The subsequent evolution, using a numerical evaluation 
of the DTQW Q is shown in Fig. 1(a), using a reason- 
ably small value of cos 9. The solution of the continuous- 
time limit p5|) can be found analytically in terms of the 
regular Bessel function ■ 



,t) ) 2 \ J„(27i)-f iJ„+i(27<) 



whose probability density is shown in Fig. 1(b). Note 
that the initial condition (|16|l differs from most previous 
studies of the DTQW: (i) it has support on n = and 
n = ±1 and (ii) it is entangled. Both are responsible 
for the slightly blurred appearance in comparison to the 
CTQW i'{n,t) = e-2^T*i"J„(27i) shown in Fig. 1(c). 
Nevertheless, the visual agreement between all three is 
quite clear, and the analytical result (|17|) shows that the 
approximation found previously |ll| is exact in the 6 
7r/2 limit. 



This result— the fact that the DTQW limits to the 
CTQW — agrees with, from an entirely new perspective, 
the interpretation of the quantum walk as a simple in- 
terference process in which the coin degree of free- 
dom is, at least in this simple case, irrelevant to the 
speedup found in quantum walk algorithms. In this 
context, it has been shown that a discrete-time evolu- 
tion can be constructed without a coin |l2l |. Here, one 
takes the Hamiltonian that generates the CTQW and 
splits it into "even" and "odd" terms H = F™'^" + H°'^'^ 
|25)]. Such a splitting is favorable for computation (clas- 



odd 



are 



sical or quantum 0), since H'^'"'^" and H 
both block-diagonal matrices, with each block a simple 
2 by 2 matrix. Then, one constructs the unitary op- 
erator U{ei,92) ^ exp(i02^f°'''')exe(«6liiJ^"'=") (the case 
9i ~ 02 — tt/A was considered in 12]). Though this oper- 
ator trivially becomes equivalent (in the continuous-time 
limit) to the CTQW for 9i = 02 ^ 0, the hmit found 
here is quite distinct. By using the even-odd splitting of 
the lattice, one can formally reintroduce the coin degree 
of freedom and recover the DTQW © with 61= 9-tt/2 
and 02 = Ti 12. It is this formulation that is relevant here, 
and the corresponding limit — > is just as puzzling. 

At this point, a number of questions arise, such as: 
what about the classical case? Is there some parallel to 
classical random walks? To answer this, it is necessary 
to identify the classical discrete-time process analogous 
to (jSJ. This is the persistent random walk |l8| ] 



PF(.{n, T +1) = apR{n - 1, r) + (3pL{n - 1, r) 
Pl(71,t + 1) = apL{n + l,r) + (3pji{n + 1,t), 



(18) 



with a + f3 — 1, where a = cos^O. This classical pro- 
cess could arise by measuring the coin after each step of 
the DTQW [11 ll9j. In general, the persistent random 
walk differs considerably from the simple diffusion pro- 
cess quantized by For example, it is known that a 
continuum limit of l|18|) with a ^ 1 {6 Q) 1261 yields 
the telegrapher's equation for p = pn + pl |l8lll9l ] 



dtpix,t) = D[dM^,t) - v-'dtp{x,t)] 



(19) 



This evolution has characteristics of both wave propaga- 
tion and diffusion, becoming the latter only in the limit 
D — > cxd; the relationship with the Dirac equation was 
explored some time ago ^] . 

To complete the correspondence with and J^J) re- 
quires consideration of the limit a — > (f? — + 7r/2), which 
even classically appears quite unfamiliar. Nevertheless, 
one can proceed by iterating H18|l once to find 



PR{n, T + 2) = a^PRin - 2, t) + af3[pL{n - 2, t) + pL{n, r)] -t- P'^PR{n, r) 
PL{n,T + 2) = a^pLin + 2,T) + aP[pR{n + 2,t) +pR{n,T)] + /3^pL{n,T). 



(20) 
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Then, letting t = rAi, a ^ 1 - (3 = 2-fAt, the hmit 
At — > leads to the result 



dtPR{n,t) 
dtPL{n,t) 



■l[PR{n - 



-2,t) 
2,t) 



Finally, defining i) = pf{{n,t)+pL{n- 
by direct subsitution that 



PR{'n,t)]. 

(21) 

l,t), one finds 



dtp{n, t) = 7[p(n + 1, t) - 2p(n, t) + p{n ~ 1, i)]. (22) 

Thus, the continuous-time limit, with a — > 0, of l|18|) 
leads to the discretized diffusion equation H22|) . in com- 
plete parallel to the quantum case. 

One might also ask: what about higher-dimensional 
walks? Is this limit restricted to one dimension? In- 
deed, beyond one dimension, or with general graphs, 
things might not be so simple. Consider, for example, 
the unitary cellular automaton proposed by Bialynicki- 
Birula [i^ . Here a four-component amplitude on a three- 
dimensional lattice '^(nx,ny,nz,T) has an update rule 
constructed from conditional displacements and coin ro- 
tations just like the DTQW. Its momentum-space prop- 
agation matrix looks quite similar to (|7|): 



U = e 



--iO<7^®I 



(23) 



but, due to the fact that operators do not commute, a 
calculation similar to JSJ immediately breaks down. This 
process depends not only on the directions a particle 
moves on the lattice, but on the actual sequence of its 
steps. A similar effect was seen in a recent quantum 
walk search algorithm ^22] (in two dimensions), where 
different coin implementations (step sequences) yielded 
different results. 

However, by using the slightly more symmetrical se- 
quence of operators 



U ■ 



xe 



— iky <7 z 



y/2p ik^ 



(24) 



the 6 ^ 71/2 limit proceeds as above ^27] leading to 

= ^4 =±1} + d^.ny + dy,nz + d^,t), 

(25) 

an obvious generalization of to three dimensions. 
Thus, the limit found here does generalize to higher- 
dimensional quantum walks. 

In summary, I have found the precise limiting proce- 
dure needed to map the discrete-time quantum walk to 
the continuous- time quantum walk. This procedure was 
extended to the classical persistent random walk and dif- 
fusion on the lattice, and to higher-dimensional quantum 
walks. In all cases one finds the initially counterintu- 
itive result that a process in which a particle moves left 
and right, flipping its direction with nearly unit probabil- 
ity per time-step, in the limit of continuous time, yields 



genuine dynamics. Classically, one finds diffusion, while 
quantum mechanically, one finds wave propagation on 
the lattice. That this connection remained hidden for so 
long, and the difficulties encountered in higher dimen- 
sions, indicate that there remains a great deal to learn 
about these simple quantum algorithms. 

I sincerely thank A. J. Dragt, P. R. Johnson, and S. 
Lloyd for helpful discussions and comments. 
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The precise limit requires introducing space-time lattice 
spacings of Aa; and At. Then, one lets 9 = mc^At/h, 
Ax = cAt, taking the limit At with the continuum 



5 



coordinates x = nAx,t = rAt held constant. 
[24] The condition ^ - (n, t) = need only be enforced at time 

t — 0. Linearity of the differential equation 1121 ensures 

that it will then be satisfied for all t 
[25] For example, with = 25„,m - (5„+i,m — (5n-i,m let 

H'^^m ~ 5n,m— + ( — l)")'5n + l,m— 1(1 — ( — l)")5n-l,m 

and H°^% = - i(l - (-l)")<5„+i,„ - i(l + 

( — l)")(5„_l,m. 

[26] The precise limit requires introducing space-time lattice 
spacings of Aa; and At. Then, one lets l3 = 1 — a = 
v^At/{2D), Ax = vAt, taking the limit At with the 



continuum coordinates x = nAx, t = rAt held constant. 
[27] Specifically, one finds that U{t) 

exp(— i$) exp(— i_fft) where <1> = tii/2 and 
H = ~~2'y{aax ® I + b ■ ay ® a), with a — 
cos^ kx cos^ ky cos^ kz, bx = sin kx cos kx cos^ ky cos^ kz , 
6y = cos kx sin fcy cos ky cos^ fcz , and bz — 
cos fca; cos ky sin ^2 cos ^2 . The two distinct eigenvalues of 
H are ±27 + &| + ^ 52 _ -|-2-y cos fcj, cos fc^ cos kz . 
Restricting the evolution to one eigenvalue, and Fourier 
transforming back yields 1251 . 



